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Abstract 

This paper is concerned with the fractionalized diffusion equations governing the law of the 
fractional Brownian motion We obtain solutions of these equations which are probability 

laws extending that of Our analysis is based on McBride fractional operators generalizing 

the hyper-Bessel operators I and converting their fractional power 1“ into Erdelyi-Kober fractional 
integrals. We study also probabilistic properties of the r.v.’s whose distributions satisfy space-time 
fractional equations involving Caputo and Riesz fractional derivatives. Some results emerging from 
the analysis of fractional equations with time-varying coefficients have the form of distributions of 
time-changed r.v.’s. 

Keywords: Fractional Brownian motion; Grey Brownian motion; Fractional derivatives; Generalized 
Mittag-Leffler functions; Riesz fractional derivatives. 


1 Introduction 


In this paper we consider fractional diffusion equations with time-varying coefficients. This investiga¬ 
tion is inspired by the fractional extension of the diffusion equation governing the law of the fractional 
Brownian motion t > 0, 0 < H < 1. Fractional Brownian motion is a zero-mean Gaussian 

process with covariance function 

Cov(B„(s),BH(t)) = -||f|^ + |sP"-|t-sp"|, t,s>0. (1.1) 

We first deal with equations of the form 

r d 

( = u„(x,t), ae(0,l), He(0,l), xeR, (1.2) 


where the operator appearing in (I1.2D is a special case of the fractional power L“ of the hyper-Bessel- 
type operator 

d d d 


T — /-“l _/-“2_ _ t^n+1 

dt dt"-dt ’ 


t > 0, 


(1.3) 


where u^, ; = 1,..., n + 1, are real numbers and n e N. 
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For a = 1, (11.21) coincides with the equation governing the one-dimensional probability law of 
while for H = 1/2, (11.211 is the classical time-fractional diffusion equation. There is a consider¬ 
ably large literature about fractional diffusion equations with constant coefficients, dating back to the 
second half of the Eighties for a detailed review see [l2lll ~). Amongst the more recent papers 

we refer, for example, to 1113.1^ [^] . In some specific cases, the relationship between the iterated 
Brownian motion and fractional diffusion equations was established and fully investigated IllL lBOf] . 

The theory of powers of Bessel-type operators was developed by Dimovski [0], McBride [l23l. [2^. 

a McBride and Lamb The fractional powers of Bessel operators were recently considered in 
in the analysis of random motions at finite velocity. 

The theory developed by McBride permits us to express the operator L“ in terms of products of 
Erdely-Kober fractional integrals defined as 




m 


r(a) 


^ t 

(t"‘ - u'"“ du, a > 0, t > 0, m > 0. 


(1.4) 


As far as we are aware, the generalization of time-fractional equations with time-varying coeffi¬ 
cients by means of the McBride theory have not been studied so far, neither from an analytical nor 
from a probabilistic point of view. 

The fundamental solution of (11.21) can be written as 


1 

t) = 2l-a/2(.Ha^-“/2,l-a/2 



X e R, 


(1.5) 


where 

00 

is the Wright function. 

Formula (II.5D shows that the solution u„(x, t) of (11.21) coincides with the fundamental solution of 
the time-fractional diffusion equation with the time-scale change t We prove that u^ix, t) is the 

law of the r.v. X„ ^(t) which is connected to the fractional Brownian motion by means of the relation 


^a.H(0=^2a,H(|BH(0|-), 0 < a < ^. 


We note that the r.v. ^(t) has variance 


^2Ha 


VarX2aH(f) =-^;-T 

’ 2“-ir(a-|-l) 

and thus the presence of the fractional derivative has a decreasing effect on the variance. 
In Section[3^we establish a relationship between solutions of fractional equations 

. a d^u 

tl-2H — 

dtj dx^ 

and solutions of higher order diffusion equations 


, , d^u 

— = (-l)''H—r, 
dt dx^ 


k>2. 


(1.7) 


( 1 . 8 ) 


(1.9) 


( 1 . 10 ) 


This type of relationship frequently appeared in the recent literature in the case of fractional and 
higher-order diffusion equations with constant coefficients (see e.g. [U). 

In the second part of this paper, we consider fractional equations of the form 


^Dl+u[x,t)=Ht^^ i|^u(x,t), ve(0,l), xeR, t>0, 

u(x,0) = 5(x). 


( 1 . 11 ) 
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where is the v-th order Caputo time-fractional derivative. This problem was partially considered 
in [13il l. The idea is to compare the two approaches (II.2D and (II.IIP . 

The fundamental solution of (II.IIP reads 


r +(X) 


U^X, t)= — 
Zn 


’ V ’ V 


( 1 . 12 ) 


where the function 


p r 1 4- V it ri ^ +0 + 1 ) 


k=l j=0 


r (a (jm -I-1 -I- 1) -I-1) ’ 


a,mA < 


such that a,m > 0 and a(jm + 1) ^ Z , 


(1.13) 


is the generalized three-index Mitta g-L effler function first introduced by Kilbas and Saigo in [j^, 15, 
M ; see also the recent monograph [Il2fl . 

Our research is related to the recent papers [HIS], where fractional generalizations of the heat 
equation for the fractional Brownian motion were inspired by some physical problems. In particular, 
in these papers the authors study 


C r,v 


Dl+u[x,t) = Dfjt 




V G (0, 1), X G R, t > 0, 


(1.14) 


in order to take into account the joint action of trapping and correlated fluctuations for anomalous dif¬ 
fusions in heterogeneous media. In this framework we recover some of the analytical results discussed 
in i]. For the special case H = 1/4 and v = 3/4, we are able to give an explicit and interesting form 
for the fundamental solution of ( I1.14P as the law of a time-changed Brownian motion. 


2 Preliminaries 


In this section we recall the main definitions and properties of fractional derivatives and fractional 
powers of hyper-Bessel-type differential operators for the convenience of the reader. There is a large 
literature about fractional powers of differential operators. We refer, for example, to the monographs 
[|23] and lUJ. We start our short review recalling the definitions of Riemann-Liouville fractional 
integrals and derivatives. 

Definition 2.1 (Riemann-Liouville integral). Let / g where —oo <a<t<b<oo, be a 

locally integrable real-valued function. The Riemann-Liouville integral is defined as 


r t 


4V(t) = 


r(a) 


/(u) 


(t-u) 


—du = if*Kjit), a>0, 


( 2 . 1 ) 


where K^it) = t“ ^/r(a). 

Definition 2.2 (Riemann-Liouville derivative). Let f g L^[a, b], —oo <a<t<b<oo, and f *K^_a e 
[a, b], m = \ a], a > 0, where [a, b] is the Sobolev space defined as 

( d™ ) 

[a, b] = |/ e Li [a, b ]: —/ G L^ [a, b] |. (2.2) 


The Riemann-Liouville derivative of order a> 0 is defined as 




(t-s)” 


7(5)d5. 


(2.3) 
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It is a simple matter to show that the Riemann-Liouville fractional derivative is rigorously the 
fractional power of the operator D = d/dt in the classical theory of fractional powers of operators (see 
e.g. [@] and the references therein). We now recall that the so called Caputo derivative, widely used 
in applications, is in fact a regularization of the Riemann-Liouville derivative. For n e N, we denote 
by AC" [a, b] the space of real-valued functions / (t), t > 0, which have continuous derivatives up 
to order n — 1 on [a, b] such that (t) belongs to the space of absolutely continuous functions 
AC [a, b] : 

AC" [a, b] = |/ : [a, b] ^ R: U) e AC [a, b] |. (2.4) 

Definition 2.3 (Caputo derivative). Let a > 0, m = \a], and f e AC'"[a, b]. The Caputo derivative of 
order a> 0 is defined as 


C 


D:.fit)=r-^—m= 


r(m — a) 


(t-s)'"-i-“-/(s)ds. 


(2.5) 


Note that in the space of functions belonging to AC™ [a, b] the relation between Riemann-Liouville 
and Caputo derivatives is given by the following 

Theorem 2.1. For f e AC'"[a, b], m= \a], a e R^\N, the Caputo fractional derivative of order a of f 
exists almost everywhere and it can be written as 


f 


^^Dl,f{t) = Dl, fW-Y, 


V 


k=0 


it-ay 
k\ 




( 2 . 6 ) 


For the proof of Theorem 2.1 consult [Q], page 50, Theorem 3.1. 

We observe that the difference between the two definitions is given by the Riemann-Liouville derivative 
of the Taylor polynomial of order m — 1 centered in the lower extremum of integration t = a. Hereafter 
without loss of generality we will consider a = 0. 

In view of this well-known relation between the two different fractional derivatives, we now recall 
the definition of fractional hyper-Bessel-t)^e differential operators, studied in a series of works by 
McBride [j^, 2£, 25] and McBride and Lamb [H^ . In [25] McBride considered the generalized h 5 ^er- 
Bessel operator 

d d 

L = t“‘ —t"2...t"" —t“"+i, t>0, (2.7) 


dt 


where n is an integer number and Uj,. 


• •) '^n+l 


dt 

are real numbers. 
The operator L defined in ( I2.7D acts on the functional space 

Fp,^ = {f:t->^fit)&F^}, 

where 


dt^ 


F = f/ eC“: t*^^ gLP(0,oo), fc = 0,l,. 


( 2 . 8 ) 


(2.9) 


for 1 < p < 00 and for any real number p (see for details [123112411 ). The operator L was first introduced 
and studied, as far as we know, by Dimovski |0] • 

Let us introduce the following constants related to the general operator L. 


n-M 




m = \a — n\, 




k = 


k=l 


P=k+1 


The definition of the fractional hyper-Bessel-type operator is given by 
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Definition 2.4. Let m = n — a>0, ^ and 


bk := i T] gK: mi 7 + M+m 7 ^ - - mZ, Z = 0,1,2,... L k = l,...,n 


Then 




k=l 


where, for a > 0 and mr\ + /i + m > ^ 


( 2 . 10 ) 


lyf = 


^—mrj—ma 

r(a) 




( 2 . 11 ) 


and for a < 0 


lyf = (tj + a + i)/;:2’“+v + f f-^/V 

m w J m j ^ m y df^ J 


( 2 . 12 ) 


For a full discussion about this approach to fractional operators we refer to [l25n . page 525. Accord¬ 
ing to the analysis developed in this paper, the domain of the operators L“ is the space F^ 

A key-role in the next sections is played by the regularized Caputo-like counterpart of the operator 
(12.101) . In analogy with Theorem l2.1l we introduce the following 

Definition 2.5. Let a be a positive real number, m = n — a>Q, f& Fp^^ is such that 

f 


L“ 


b-l 




k=0 


exists. Then we define ^ L‘^ by 


^L'^fit) = L'^ 


b-l pJc 




k=0 


(2.13) 


where b = [a]. 

The relevance of this definition for the applications is due to the fact that for physical reasons we 
are interested in solving fractional Cauchy problems involving initial conditions on the functions (and 
their integer order derivatives), while in the case of Riemann-Liouville operators we must consider 
initial conditions involving fractional integrals and derivatives. 


3 Fractional diffusions with time-dependent coefficients 

Let us consider the fractional diffusion equation with time-dependent coefficients 

^^JtJ ^a(.x,t)=H'^—u„ix,t), (3.1) 

where H e (0,1) is the Hurst parameter, a e (0,1) and, according to Definition 2.4, the fractional 
operator is defined (for n = 1, Uj = 1 — 2H, a 2 = 0 and thus bj = 0, m = 2H) as 

= <^3.2) 
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where “ is the Erdelyi-Kober fractional integral defined in (12.IIP . 

According to lUJ, the fractional power of the operator admits the representation (I3.2P 

in the domain fp where p e R should be taken such that 


= 0 C „ 2h : e K : 2 Ht] + fX + 2H^-- 2HI, 1 = 0,1,...}. 

p 


This constraint gives the implicit relation between the parameter p of the domain of the operator and 
the parameter H of the specific operator that we are considering here. 

We can call (13.IP a bifractional equation, depending on the Hurst parameter H (which, in the case 
a = 1, is related to the fractional Brownian motion) and the real parameter a (that is the fractional 
power of the operator). Notice that for a = 1 we retrieve the heat equation with one time-dependent 
coefficient. 


du d^u 
=H- 


dt 


dx^ ’ 


(3.3) 


satisfied by the probability law of the fractional Brownian motion. 

In view of definition 2.5 the regularized Caputo-like counterpart of the fractional operator (13.2P 
reads 

r r t“^“ 

u„(x,t)=(^ti-2»—J u„(x,t)-(2H)° ^^^_^^ u„(x,0), (3.4) 

in complete analogy with the classical theory of fractional derivatives. 


Theorem 3.1. The solution to the Cauchy problem 


|C aG(0,l), XGR t>0, 

|u„(x,0) = 5(x), 


is given by 


1 f 2^l'^\x\ 

Ua{x, t) = 2i-a/2(.Ha'^-“/2.1-a/2 I pHa 


Proof. By using 03.4P we obtain the equation 


A-2H 


dt 


' t 

u„(x, t) = H“—^u„(x, t) -I- (2H) 
dx^ 


r(l - a) 


5(x). 


The Fourier transform of 03.7P reads 


) ii„(/3,t) = -H“^2ii„(/3,t)-|-(2H)“ 


i.-2Ha 


r(l-a) 


The solution to 03.8P can be easily determined as 


u^iP, t) = 


/32f2 


2 “ 


where Ea_i(0 is the Mittag-Leffler function. Hence, by using 03.2P . we have that 
3 


A-2H 


dt 


ii„(/3,t) = (2H)“t 


-2HajO,-a 

■'2H 


___ , a2f.2Ha 

k=0 


T 00 / a2\^ f2Hak-2Ha 

(™)“Z -k 


r(afc-bl) 


k=0 


(3.6) 


(3.7) 


(3.8) 


(3.9) 


(3.10) 


r(afc -I-1 — a) 


6 













\ t 


r(l - a) 


= -H“/ 3 ^u„(^,f) + ( 2 H)“ 


r(l-a)’ 


where we used the fact that 


jQ-a 2aHk _ r(afc + 1) 2aHk 

r(afc + l-a) ■ 


By inverting the Fourier transform u^P, t) we obtain the claimed result. 


(3.11) 

□ 


We note that the time-change t reduces equation (13.11) to the classical time-fractional diffu¬ 

sion equation 

1 

■ ' (3.12) 


with fundamental solution given by (see e.g. fl3^ ) 


Ua(.X,t) = 


2l-a/2(-f 


i'^-a/2,l-a/2 




t 2 


(3.13) 


This is a particularly interesting fact. In the first part of the paper in fact we prove formally that 
suitably time-changed fractional diffusion processes can be analyzed by means of the McBride theory. 
The converse is not true. More general equations can be treated by applying the same theory but 
cannot be recovered by means of simple deterministic time-changes. 


Theorem 3.2. The solution to 

fc ae(0,l), xeK, t > 0, 

lu„(x,0) = 5(x), 


can be represented as 


Ua{.X,t) = 


hV2 

sprit^ 






’^2aix,z)dz, 


where is the solution to 


and 


^ “^2aix,t)=H^^f^U2a<ix,t), a £ (O, i] , 

P 2 a (^.0) = 5(x), 

(l^“^"|^)^“u2„(x,t)=H2“|iu2„(x,t), ae(i,l), 
lt 2 a(^> 0 ) = 5(x), 
yj-p2aix,t)l^O = 0. 


(3.14) 


(3.15) 


(3.16) 


(3.17) 


Proof. The proof follows the main steps of Theorem 2.1 in [13011 . In more detail, by applying the 
duplication formula of the Gamma function we have 


r(-—-Fl--) = /7r2“^''+i^ 


2 2 

By substituting (13.181) in (13.61) we have 


,r(l-a(fc-Fl)) 

I ^ p _ a()c+l) ^ 


(3.18) 


^aix,t)= 


k=0 


fc!/7r2“^‘^+^^“^r(l - a(fc-F 1)) 


(3.19) 
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2l+a/2 

1 

2l+a/2 


2l+a/2y^ff: 


E 

k=0 


(“7^) 

fc!2^r(l-a(fc + l)) 


-a/2-1/2 


(-|x|/(wt2«)“/2)*' 


S 


dw 




1 r“ 1 

=—= e~'^w~^^^U2a{x,^ft(2^w)'^)dw 

V^Jo 

= [z = \/7(2^w)iH ] 


Ha/2 

-/nt^ 




"5('^) i^. 


2a(x,z)dz. 


Remark 3.1. In the special case cl= h the solution reads 


□ 


HV2 , 


1 /• 

■i(^) u. 


(x,z)dz, 


(3.20) 


where Ui(x, t) is the fundamental solution of the heat equation 


t ^ )ui(x,t)=H^^Ui(x, t). 


(3.21) 


Therefore 0 coincides with the probability density of the r.v. Bjj where Bjj and B^ 

are independent fractional Brownian motions. We note that YarB^^t) = t^^ and YarB^(t') = 4f^^. 

If H = ^, we recover the result discussed in ^3^1 . i.e. Ui/ 2 (.x, t) coincides with the probability density 
of the iterated Brownian motion B^ QB^(t)Q (YarB^{t)= t andYarB^{t) = dt) . 

Remark 3.2. We observe that, in the more general case, for a e (0,1)^ the solution u^Cx, t) coincides with 
the probability density of the r.v. X 2 a,H where X 2 a,H^^')> t > 0, is a r.v., independent from 

Bu^t) with density law equal to U 2 a{.x, t). 

Remark 3.3. It is simple to find the moments of the process Xg^u^t), t>0, by means of its characteristic 
function ( 13:^ . In particular, its variance is given by 


VarX„,„(t) 


^2Ha 

2“-ir(a + l)‘ 


(3.22) 


We observe that, for a = 1 in (I3.22D . we recover the variance of the fractional Brownian motion; for 
H = 1/2 we have the variance of the time-fractional diffusion process. 

More generally, from (IT^ . we infer that 


EXf„(t) = 


r(2m + l) 

2™T(arn + l)‘ 


m e N. 


(3.23) 


Remark 3.4. 

coefficients 


We now consider the more general time-fractional diffusion equation with time-varying 


u„(x,0 = 


1-a^Y 

—^U^[x,t), Gj < 1. 


(3.24) 
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In view of the McBride theory we have that in this case 




The solution of the Cauchy problem 

1 '^ 8 ^^<^<ix,t), aG( 0 ,l), 

|ii„(x,0) = 5(x), 


is given by 


Ua(x,t) = 


2l-a/2(-—a 


i-a, ^-a/2,l-a/2 




Remark 3.5. If we consider the complete McBride operator of order one, thatis [t‘‘^-^t'^^), witha-j 
1, it is easy to prove that the function 


h„(x, t) = 


solves the non-homogeneous pde 


-W_ 


. ... l-ai-a 2 ” 0 !/ 2 ,l a/2 1 1 -dj-ci 2 _ 

2l-a/2(-—5—a+aj ^—5-c 




dt 


fOi (.aj y^^(^y.,t) = ^^hJ^X,t) + 


3x2 


5(x) 

(1-01-02)““ r(l-a)' 


Indeed, by taking the Fourier transform of ( I3.29D we have that 


KW,t) = -pXw,t) + 


^-a[l-ai-a2)-a2 ^ 


whose solution is given by 


KW,t)=—E^_, ( - 


(I-O1-O2) “r(l-a)’ 

p2 i.ail-a^-a2) 


(1 - Oi - 02 )“ 


In view of the fact that 


ffj- I ^ 1 AJ 1 -a,-a. 


= ( 1 - 01 - 02 ) 


a y l-ai-a2y 


^r-aCl-Qi-aj) 

*- J 


we have that 


f d ^ 

\“ 1 ( 

I i-^1 _ f<^2 

I 7sr^a,i ( 

V dt , 


p2j.a[l-ai-a2} 
(1 - Oi - 02 )“ 


( 1 - 01 - 02 )“^ 
00 

-s 


(-P") 


2 'vA: f.a{l—ai—a2)k—a[l—ai—a2)—(i2 


(1 - Oi - 02)“^^ r(afc + l-a) 

^_-^'^k^P'^2k+2 ^aCl-aj-Oj/Jc-aj 


(1 - Oi - 02)“^^ r (afc + 1) 

j-a(l-ai-a2)-a2 ^ 


--P “a + 


(1-01-02) “ r(-a + l) 


(3.25) 

(3.26) 

(3.27) 

+ 02 < 

(3.28) 

(3.29) 

(3.30) 

(3.31) 

(3.32) 

(3.33) 
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By inverting the Fourier transform we obtain the claimed result. The case a-^ + 02 = 1 can be treated by 
using the results discussed by Lamb and McBride in Hj. In this case, a logarithmic deterministic time- 
change leads to the solution. Finally, in the case aj + 02 > 1 complex coefficients appear and therefore we 
neglect this case. 

If we consider the space-time fractional equation 

ve(0,2), xeR (3.34) 

where 3''ld\x\'' denotes the Riesz fractional derivative 


QV 


gU) 


1 15™ 

r(m — v) 2 cos 



\x — s 




T^g(5)d5, 


m — 1 < V < m, X > 0 (3.35) 


where m is the ceiling of v, we have that its fundamental solution reads 


ii„(x, t) 


1 

sin a 71 


Ea,l - 


\P\M 


v,.2Ha 


2 “ 

..a—1 


dp 


' 2r“ cos an -I-1 


g-i^Xg-|^;r/“V'-d^dr. 


(3.36) 


The expression u„(x, t) coincides with the distribution of the time-changed r.v. Q#'„(t)^ where 
is a symmetric stable r.v. of order 0 < v/a < 2 and #'a(t) is the Lamperti r.v. with parameter 
a (that is the ratio of two independent positively skewed stable r.v.’s of order 0 < a < 1 ), possessing 
density 

sinna 


Pi^a e dr}/dr 


71 r^“-f 1-f 2r“ COS Tia’ 


r > 0, a e (0,1). 


(3.37) 


3.1 Relation with the generalized grey Brownian motion (ggBm) 


The generalized grey Brownian motion (ggBm) was recently introduced and studied by A. Mura and 


2811291] as a family of non-Markovian stochastic processes for anomalous fast or slow 


coauthors in 
diffusions. 

The marginal density function of the ggBm is given by 


P(x, t) = 


-i-£ 

2tT/2 2’-^ 2 


_ m 

tr/2 J ’ 


re(0,2], 5e(0,l]. 


(3.38) 


The ggBm includes for 5 = 1 and y e (0,2] the fractional Brownian motion, for 0 < 5 = y < 1 the 
grey Brownian motion introduced in the literature by Schneider OSfl and for 7 = 5 = 1 the standard 
Brownian motion. We observe that the density (I3.38D coincides with (13.61) for 5 = a and y/2 = Ha (up 
to a constant due to the scale invariance of the Wright function). This means that the equation (13.51) 
governs a ggBm of parameters (a, Ha) explicitely dependent on the Hurst parameter of the fractional 
Brownian motion. 

In his PhD thesis [|26|] . A. Mura introduced the fractional equation governing the probability density 
of the ggBm, that is given by 


P(x,t) = P(x,0)-l- 




r(5)5 


(3.39) 


Equation ( I3.39D is an integro-differential equation involving the Erdelyi-Kober integral. 
We now show the equivalence between the fractional equation 

3 3^ 

3t I 3x^ 


10 


(3.40) 



















and the master equation (I3.39I) introduced by A. Mura, in the case of 5 = a and y/2 = Ha. 

From (12.131) and by using (13.21) . we have that 

P(x,t)= f ti-2«—J (P(x,t)-P(x,0)) (3.41) 

= (2H)“t-2«“7^-“ [P[x, t) - P[x, 0)) 

3^ 

= H“^P(x,t), 

3x^ 


and therefore, 


4 “(P(x,t)-P(x,0))= 


.2Ha 


2“ 3x^ 


P(x, t). 


(3.42) 


We now recall the following property of the Erdel 5 d-Kober integral (see 1 I 25 I 1 . Theorem 2.7, page 523) 




(3.43) 


Thus, in our case, we have that the inverse of the operator appearing in the left hand side of ( I3.42D is 
given by 

(3-44) 

Finally, by applying the inverse operator ( I3.44D to both sides of ( I3.42D . we arrive at 


P(x,t)-P(x,0) = 


-‘ 2 H 
2 “ 

2i-“H 

r(a) J 


.2Ha 




(3.45) 


32 




which coincides with (I3.39D for 5 = a and j = 2Ha, up to a multiplicative constant. 

A further approach highlighting the relation between ( I3.39D and equations involving Erdelyi-Kober 
derivatives was recently discussed by Pagnini in [l33l l . To conclude, we remark that the general McBride 
theory for fractional powers of hyper-Bessel-type operators represents a convenient framework to de¬ 
rive the governing equation of the ggBm considered as a fractional generalization of the diffusion 
equation governing the fractional Brownian motion. 


3.2 Relation with higher order heat equations with time-varying coefficients 

In this section we consider the relation between the solutions of Theorem 3.1 and the fundamental 
solutions of higher order heat-type equations. 

We start from a special case of (I3.5D . corresponding to a = 2/3. We know, from [|3^ (Remark 4.1, 
pages 231-232), that the fundamental solution to the fractional equation involving Caputo derivatives 




3x2 ’ 


A > 0, t > 0, 


(3.46) 


is related to the fundamental solution of the linearized Korteweg-DeVries equation, that is the third 
order heat equation 

3v „d^v 

— = -A2—, A>0, t>0. (3.47) 

3 1 3x'^ 

By following the same idea, we establish here a relationship between fractional equations and higher- 
order heat-type equations with time-varying coefficients. 
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Theorem 3.3. The solution U 2 /k[x, t) of (I3.5D . for any fc e N ; fc > 2, coincides with the fundamental 
solution v(x, t) of the higher order diffusion equation with time-dependent coefficients 


dv , d^v 


X > 0, t > 0. 


dt '' ^ dx^ 

Proof We first recall that in the special case a = 2/k the solution of (13.51) is given for all x e R by 

1 (-1)" 2"/‘^|x|" 


(3.48) 


Ui/kix, t)= — 

O 


■s- 


Hn 


2i-u2H/;c^n!r(-| + l-l) ff 
By direct calculation we have, for x > 0, that 

^11 f ']^ri»~)n/k^ni. — THn—^H—2H 


dt 


n=0 




and this coincides with (—1)*'H times the fc-th order derivative with respect to x, as claimed. 


(3.49) 


(3.50) 


□ 


Remark 3.6. As a special case of interest of Theorem \3.3\ we observe that the solution U 2 /-^{x, t) of (13.51) 
coincides for x > 0 with the fundamental solution v(x, t) of the linearized Korteweg-DeVries equation 
with a time-dependent coefficient 

dv d^v 

t^-^— = -H^. (3.51) 

d t dx"' 

In particular, it is possible to prove that ( I3.49D can be represented in terms of Airy functions, following the 
same arguments of Theorem 4.1 in ^3^1 . since the fundamental solution to (13.511) can be written as 


f 


v(x, t) = 




-Ai 


2H 


•f¥“ 


(3.52) 


Observe that the transformations t' — t^^ and x' = x V2, reduce equation (I3.51D to 


dv 


d^v 


whose fundamental solution is 


and thus we have (I3.52D . 


, , 1 I X' 

v(x , t ) = Ai 


V 


4 Diffusion equation with time-dependent coefficients involving 
Caputo derivatives 

In this section we adopt a different fractional generalization of the time-var 3 dng diffusion equation, 
involving Caputo time-fractional derivatives. Our aim is then to compare analytical and probabilistic 
results relative to these two different approaches. 

In more detail, here we consider the time-fractional diffusion equation 


('^Dl+u[x,t)=Ht^" ^^u(x,t), vg(0 ,1], xeM, t>0, 

|u(x,0) = 5 (x), 

with 0 < H < 1. The derivative operator appearing in ( 14.ID is the Caputo fractional derivative. 
For the utility of the reader, we first recall the following result (see E3] , pag.232-233). 


(4.1) 
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Lemma 4.1. The function 


k-l 


fit) = 1+1,1 =1 + 2 (Ar+r)'^ ]“[ 


r(v; + rj + r + i) 


k=l 


LJ rCvj + yj + r + v + i)’ 


(4.2) 


is a solution to the Cauchy problem 


i^D^Jit) = At7(t), t > 0, V e (0,1], r > -V, 

1 /( 0 ) = 1 . 


(4.3) 


Proof We have that 


k-l 


k=l j=0 


r (vj + rj + r +1) 

=0 r (vj + rj + r + V +1) 


(4.4) 


^_i2^v + fcr + i)_^fc(v+r)-r]-[ r(vj + jr + r + i) 


= Va‘ 

^ r (fcv + fcf +1 - v)" j J r (vj + rj + r + V +1) 

= At’' y A''t''(''+’') ^"(''(v + r) +v + r + i) A r(vj+jr + r + i) 

^ r(r(v + r) + r + i) jjr(vj + rj + r + v + i) 


= At’' 


1 + Vifi r(vj + rj + r + i) 

^ r (vj + rj + r + V + 1 ) 


= At’'/(t). 


□ 


The function 04.2D is the generalized Mittag-Leffler function first introduced in [14, 15, 1^. An 
application of this result in fractional elasticity has been recently discussed in [ 0 ] • 

We are now ready to state the following 


Theorem 4.1. A solution u(x, t) to 04.ID can be written as 


uix,t)=— e ‘^^£:^,i+ 2 i^, 2 i^ ^)d^. 

^ J —00 

Proof By taking the Fourier transform of D4.1D we obtain the differential equation 

t) = -Ht2»-i/32[/ (p,t), t > 0, V e (0,1], 

Uip,0) = l. 

Equation ( I4.6D is solved by (see Lemma R~T] and [fl^ . page 233, examples 4.11 and 4.12) 

UiP, t) = . 

Moreover, the solution is proved to be unique for 1/2 < H < 1 (see E3], page 233). By taking the 
inverse Fourier transform we obtain the claimed result. □ 

Remark 4.1. An interesting case of ( 14.7D is given by the choice H = 1/4, v = 3/4. In this case, from 
(I4.23D we have that 


(4.5) 


(4.6) 


(4.7) 


w.t)=Ei,i,_i ) = 1+i: () n 


k=l 


^2(.i/4y^r(i + i) 


22 


4 2 ^ 

y ni + 1) 


(4.8) 
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Since 


k-i pri + 1') 


r(i/2)r(3/4) 

rk+2-^T^rk+3^T'rk+4^ 


)=0 r(i +1) r(^)r(^)r(i^) 

and using the well-known multiplication formula for the Gamma function 

r(z)r j r j ... r ^ j = (27r)^"'“i^''2m2-'"^r(mz), 


we have that 


Therefore, 


jr4 +1) _ r(^)r(3/4) 
Ur(i + ^)” 712-2'^+!fc! 


u(p.t}= E.,_, f-= 1 + 


32.1/4^*= V 4 2 


= 1 + 


r(3/4) 


7tV2 
r(3/4) 


e-'^w-^/^ 


Finally, we have 


nV2 Jo 


u(x,t) = 


) k=l 


k=i 

2^-r.-dw 


k\ 


k\ 


r(3/4) 

71 a/2 


V^r(3/4) 

Ti^ft 


g 4Vtiv 

— dw 


\/4t 


0 


;e i-ttdw. 


that is the probability law of the r.v. BiW^), where 

V^r(3/4)e" 


P{Wf e dz}/dz = 




z > 0, 


and B is a standard Brownian motion. 


(4.9) 


(4.10) 


(4.11) 


(4.12) 


Remark 4.2. For v = 1 we show that the solution u(x, t) in equation (I4.5D coincides with the density 
function of the fractional Brownian motion for a fixed time t. Indeed, for v = 1, formula (1^ yields 


U[p, t) = exp 




(4.13) 


which is the characteristic function o/B^(t), as can be checked from ( I4.2D for A = —H/3^ and y = 2H — 1 
(or directly from ( I4.3D ). 

Remark 4.3. For 1 + = 0, the characteristic function ( I4.7D becomes 


U(l3,t) = 




HI3^r(l-v) < 1, He 0, 


(4.14) 
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and is independent from t. Indeed, in this case we have that 






;t=i 


;=0 


r(i) 


00 -I 

y - v))^ = -r- 


when HjS^ril — v) < 1. 

Wfe observe that, even if the function (I4.14D is independent from t, it satisfies the time-fractional 
differential equation 

involving the Riemann-Liouville derivative. We recall that the Riemann-Liouville derivative of a constant 
function does not vanish and in more detail 


const, t 

DZ.+const. = —;-- 

° r(l-v) 


(4.16) 


Remark 4.4. For v = 2H — 1 = ^ and therefore H = |, the characteristic function ( I4.7D becomes 


2-^4" 


uw. t)=1+g(-Hp»t)‘n y j=1+g n i [7 j ■ 


1 f2j 


We observe that the coefficients in ( I4.17D have the following probabilistic interpretation 


uU"!]=Hr\ 




;=i 




(4.18) 


In ( I4.18D by ©■' we denote independent binomial r.v.’s with parameters 2j and 1/2. For large values of j, 
we know that P{Q5 (^2j, = )} = Ij and thus the coefficients of ( I4.17D vanish as l/\/H. 

4.1 Higher order extension 

Here we consider the higher-order time-fractional equation 


t) = cjf2Ht^-^ f^f(x, t), V e (0,1], 0 < H < 1, 


/(^,0) = 5[x), 


(4.19) 


and k is an integer greater than 2. The coefficient cj^ = (—jf jg gygjj ci^ = ±1 otherwise. 
The fractional derivative appearing in ( I4.19D is the Caputo fractional derivative. 

Theorem 4.2. The solution f(x, t), x e K, t > 0, to ( I4.19D is 


f(.x,t)= — 
Zn 


(c;t2H(-i/3)'^t''+2«-i) dp. 


(4.20) 


Proof. Applying the Fourier transform 


FiP,t)= e'^V(x,t)dx, 


(4.21) 
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to ( I4.19D . we have 


^Dl,F il3, t) = {13, t) , 

F(/3,0) = l. 


The solution to (I4.22D is immediately obtained: 

F(/3, t) = (c,2H(-i/3)''t"+2H-i) , 


(4.22) 


(4.23) 


Recall that the solution is proved unique for 1/2 < H < 1. From (I4.23D it is immediate to obtain 

(IHQI). □ 


In the special case H = 1/4, v = 3/4 and k = 2n, the Fourier transform of the solution of ( I4.19D 
reads 


r(-)A /2 

F(^, t) = —^77^1 


Jo 


(4.24) 


Thus the inverse Fourier transform coincides with the law of a pseudoprocess of order 2n evaluated at 
the random time W^. In a similar way the odd-order case can be treated. For a wide discussion about 
pseudoprocesses we refer for example to ifl^ and the references therein. 

4.2 Space-time bifractional equations 

Theorem 14. II can be modified by introducing the Riesz fractional derivative, in the right hand side of 
04.ID . thus obtaining 

VG(0,1], 0<a<2, 0<H<1, 

\/(x,0) = 5(x). 

This equation has been object of recent interest within the framework of anomalous diffusion processes 
0,13] . Similarly to Theorem 14.21 we apply the Fourier transform to 04.ID . arriving at 


t) = -Ht 2 H-i|^|“F(^, t), V e ( 0 , 1 ], 0 < a < 2 , 
F(^,0) = 1. 


The solution to ( 14.26D is thus 


2n 


1, 2h^ 2h^ d/3. 

> V ’ V 

D 


(4.26) 


(4.27) 


We observe that ( I4.27D coincides with equation (63) of [Q] with a suitably adaptation of the notation. 

Remark 4.5. When H = 1/4, v = 3/4, on the basis of the results discussed in Remark |TT| (see equation 
04.1 ID ), we have that 


F[p,t} = 


2y2r(3/4) 

nVt 




that coincides with the probability law of the r.v. Y„( 2 irf), where 

( 1 , 2y2r(3/4)e“^ 

T®-^ *}/“" =-Tyf— 

and is a symmetric stable r.v. of order a. 


z > 0 


(4.28) 


(4.29) 
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Remark 4.6. When H = 1/2, the solution ( 14.271 ) reduces to 


f-l\P\-tAdp. 

A \ ^ J 

The expansion of the Mittag-Leffler function yields 


fix, t) = — 
zn 


(4.30) 




^-ipx 


SinVTT 


R 

sin V 71 


— 1 


e ( 2 ) ''^drd^ (4.31) 


0 

.V —1 


+ 2 r'' cos VTi + 1 


2tc^ Jo r^''+ 2r ''cos vti + 1 




Then f{x, t) is the distribution of the time-changed r.v. Yg^i„ where Wyit) is the Lamperti r.v. 

and Y^/y is a symmetric stable r.v. of order 0 < a/v < 2. 


Remark 4.7. Ifv = jr], we can expand the Mittag-Leffler function appearing in ( I4.30I I with an alternative 
formula (see IsZl pace 2911. i.e. 




sin 771 


ri-t 

- E 1 (—r'&^^'^t'^^ dr, 

+ 2r'^cos77i + 1 ^ ■' 


where 7 e (0,1], 17 e (0,1]. In our case, for H = 1/2 we obtain that 


fix,t) = 


smyn 


271 




n 

rt-l 


rt-l 


271^ Jo + 2r’'cos77r + 1 


0 Hr+ 2 ri^ cos 771+1 ^ '' ' ^ ' 

dp dr 


(4.32) 


(4.33) 


sin 771 sinTjTi 
2n^ 


r-r-i 


dr 


r^r + 2rr cos77r + 1 J_ 


+00 


-'^^d/3 


y ^ 

0 + 1 + 2 y’i cos TiTj 


dy. 


Therefore, for H = 1/2 and 0 < a /717 < 2, /(x, t) is the distribution of the time-changed random variable 
'^a/yri where Wit) is the Lamperti r.v. and Yg^/y.^^ is a symmetric stable law of order ■^. 

Remark 4.8. We observe that for v = 1, a e (0,2) and H e (0,1), from (I4.27D we obtain the law of a 
symmetric stable process with a deterministic time-change, whose characteristic function is given by 


F(P,t) = exp\ - —l^r 


(4.34) 
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